Tunneling Anisotropic Magnetoresistance of Helimagnet Tunnel Junctions 
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We theoretically investigate the angular and spin dependent transport in normal-metal/helical- 
multiferroic/ferromagnetic heterojunctions. We find a tunneling anisotropic magnetoresistance 
(TAMR) effect due to the spiral magnetic order in the tunnel junction and to an effective spin-orbit 
coupling induced by the topology of the localized magnetic moments in the multiferroic spacer. The 
predicted TAMR effect is efficiently controllable by an external electric field due to the magneto- 
electric coupling. 
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Introduction.- Transport across two ferromagnetic lay- 
ers separated by a tunnel barrier depends in general on 
the relative orientation of the layers magnetizations^, giv- 
ing rise to the tunnel magnetoresistance (TMR) effect^. 
In the presence of spin-orbit interactions TMR becomes 
spatially anisotropic^^—. Tunnel anisotropic magnetore- 
sistance TAMR is observed not only in magnetic tun- 
nel junctions with two ferromagnetic electrodes'^ but also 
in ferromagnetic/insulator/normal- metal systems such as 
Fe/GaAs/Au^. Here we show that TAMR is a distinc- 
tive feature of normal-metal/multiferroic/ferromagnetic 
heterojunctions with the particular advantage of being 
electrically controllable. The coexistence of coupled elec- 
tric and magnetic order parameters in multiferroics^ 
holds the promise of new opportunities for device 
fabrications^^. Our interest is focused on helimagnetic 
multiferroio^iiS. The topology of the local helical mag- 
netic moments in these materials induces a resonant, 
momentum-dependent spin-orbit interaction-. The non- 
collinear magnetic order together with the induced spin- 
orbit coupling result in uniaxial TAMR with b. C2v sym- 
metry. These two factors and their interplay determine 
the size and the sign of TAMR. In particular, a linear de- 
pendence on the spiral helicity results in an electricallyii 
tunable spin-orbit interaction^ by means of the magneto- 
electric coupling, and thus TAMR is electrically control- 
lable accordingly. 

Device proposition.- The proposed multiferroic tunnel 
junction is sketched in Fig[T] It consists of an ultrathin 
helical multiferroic barrier sandwiched between a nor- 
mal metallic (NM) layer and a ferromagnetic conductor 
(FM). The ferroelectric polarization P in the multifer- 
roic barrier creates in general surface charge densities 
±|P| which are screened by the induced charge at the 
two metal electrodesi^. A depolarizing field emerges in 
the barrier. Taking the spontaneous electric polariza- 
tion as Pz ~ 700 fiC/m^ and the dielectric constant to 
be e = 30 in the ferroelectric phase of TbMnOa^, the 
potential drop generated by the depolarizing field in the 
multiferroic barrier is estimated to be on the energy scale 
of meV, which is much smaller than any other relevant 
energy scale in the system. In the present study, we ne- 
glect this potential modification, and assume that the 
barrier potential has a rectangular shape with the height 



Vq. All energies are given with respect to the NM Fermi 
energy Ep . Based on this approximation, the Hamiltoni- 
ans governing the carrier dynamics in the two electrodes 
and the oxide insulator have the following form, 



Hnm — 
Hmf — 
Hfm = 



2ni. 
2m' 
2m, 



-V^ for z<0, 

+ Jrir • cr + Vo, for < z < d, 
- - Am • (T, for z > d, (1) 



where Vq and d are the height and the width of the 
potential barrier (see Fig[ljc)), is the free-electron 
mass, m* is the effective electron mass of the oxide ( 
m*/me ~ 10), and cr is the vector of Pauli matrices, 
m — [cos 0, sin 0, 0] is a unit vector defining the in-plane 
magnetization direction in the ferromagnet with respect 
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FIG. 1. (Color online) (a) Schematic diagram of the helical 
multiferroic tunnel junctions consisting of a normal metallic 
layer as the bottom electrode and a ferromagnetic layer for 
the top one. The vector m indicates the magnetization ori- 
entation specified by the angle in xy (FM) plane. The zx 
plane refers to the spiral plane of a multiferroic oxide, (b) The 
arrows show the induced resonant spin-orbit coupling, qkxO-y 
in the multiferroic barrier, (c) The tunnel barrier potential 
profile. 
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to the [100] crystallographic direction. A is the half- 
width of the Zeeman sphtting in the ferromagnetic elec- 
trode. Jrir is the exchange field, where rir is given by 
the multiferroic oxide local magnetization at each spiral 
layer (labelled by a integer number I) along the z-axis^i, 
i.e., rir — (— 1)' [sin^r, 0, cos 6*^] with 9r = Qm • r and 
— [Qi 0, 0] being the spiral spin-wave vector. The 
physical picture behind the term Hmf in eq.([T]) is that 
a tunneling electron experiences an exchange coupling 
at the sites of the localized, non-coUinear magnetic mo- 
ments within the barrier. In effect this acts on the elec- 
tron as a non-homogenous magnetic field. Performing 
a local unitary transformation within the barrier~, one 
can also view the influence of the barrier as consisting of 
two terms a homogeneous Zeeman field, and a topology- 
induced spin-orbit coupling SOC that depends solely on 
the helical magnetic ordering. As shown in^, this SOC 
depends linearly on the electron wave vector and on the 
helicity of the magnetic order— and is explicitly given by 

SOC-- qk^ay. (2) 

2m* ^ ^ ^ ' 

The dependence on kx resembles the resonant semicon- 
ductor case when the Rashba^^ and Dresselhaus^'* spin- 
orbit interactions have exactly equal strengths. Provided 
that the in-plane wave vector k|| is no-zero, an electron 
in the oxide undergoes an exchange interaction with the 
local spiral magnetic moment and the induced spin-orbit 
coupling. So an electron spinor in the multiferroic barrier 
is determined by the following spin-dependence term, 

H^fJ ^w{0r,k)-a (3) 

where 

v^f(6'r,k) = [{-ly J sin 9r,qkx,{-^y J cos 0r] (4) 

and q = i^^q- With this effective spin-orbit interaction 
we analyze the angular dependence of the electron tun- 
neling through the helical multiferroic barrier. 

Phenomenological theory.- We assume the strength 
of the effective Zeeman field |w(6'r,k)| to be relatively 
smaller than the Fermi energy Ep and the band split- 
ting A. Proceeding phenomenologically as in Refi one 
expands the transmissivity as a perturbative series of 
m • w(0r,k). Up to the second order, the transmissiv- 
ity reads, 

r(k,m) =a(°)(k) -fa^')(k)[m-w(0„k)] +afV(^r,k)|2 
-|-4''(k)|m-w(0„k)|2 (5) 

The expansion coefficients, a'p [i — 1,2; j = 

0,1,2) satisfy the symmetry relations a^^\kx,ky) = 

flj ( kx, ky), Oj {kx,ky) = ( kx,ky) and 

af\kx,ky) = af\ky,kx). Based on the linear-response 
theory, the conductance G'((/)) is found as 



G{4>) ^-f^ j m) = Go + Gamso{4>) (6) 

where Go is the angular-independent part of the conduc- 
tance, and 

Ganiso 

(0)=5tr[AM(</))] (7) 

is the anisotropic spin-orbit coupling contributions and 
g = e'^/Sn^h. A and M{<t)) are matrices whose elements 
are given respectively by 

Aij = {a\ '{k)wiWj), Mij{(f>) = mim-j = x,y,z). 

(8) 

The notation (...) stands for the integration over Or and 
k||. Introducing w(0r,k) and m into Ganiso{4>), Eq-©, 
the anisotropic conductance can be rewritten as 

Ganiso{4') = Oi COS^ 4l + 13 SIT? (j). (9) 

Considering the symmetry of the expansion coefficient 

flj (k) we obtain for the above expression, a — 

g{a^^\\i)j'^sin^ dr) and /? = g{a^^\\i.)q^kl) . Hence, the 
TAMR coefficient is given by 

TAMR = - 7(1 - cos2^), , = ^10) 

The above angular dependence of TAMR is quite general. 
The helical magnetic order and the induced spin-orbit in- 
teraction give rise to the anisotropy in the magentocon- 
ductance. However, as evident from Ea. ([TU)) . the TAMR 
coefficient 7 depends on (a — /?), i.e. contributions from 
the exchange interaction and the spin-orbit coupling term 
have opposite effects, which is confirmed by the following 
model calculations. 

Ultrathin Carriers. -Experimental observations^ indi- 
cate that thin film multiferroics can retain both magnetic 
and ferroelectric properties down to a thickness of 2 nm 
(or even less) . To get more insight in TAMR we consider 
ultrathin tunneling barriers that can be approximated by 
a Dirac-delta functiooi^. The effective spin-orbit inter- 
action H^jp throughout the multiferroic barrier reduces 
then to the plane of the barrier, H^jp — \v{9r, k) • a5{z) 
with w(0,.,k) = [J sinOr, qkx, J cos9r]. J and q are 
renormalized exchange and resonant spin-orbit coupling 
parameters, q ss qVod and J ~ {J{z))d referring to space 
and momentum averages with respect to the unperturbed 
states at the Fermi energy. In the following, we treat J 
and q as adjustable parameters. Obviously, the electron 
momentum parallel to the junction interfaces k|| is con- 
served. Then the transverse electron wave functions in 
NM (z < 0) and FM {z > 0) regions can be written as 
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with 




K = \KE + aA)/ 



and the spinors introduced as 



Xa 



1 
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(12) 



(13) 



(14) 



(15) 



and correspond to an electron spin parahel (u = 1) or an- 
tiparallel (cr = — 1) to the magnetization direction in the 
ferromagnetic electrode. The reflection (ro-^o- and Va-^a) 
and transmission {ta,a and ta,a) coefficients can be ana- 
lytically obtained from the continuity conditions for ^'(z) 
and vl/'(z)/m at z = 0^, 



nM(O-) - *?^m(0+), 



(16) 



2me 



dz 



+ (Vbd + w(0,,k).^)*^M(O") 



2me dz 



z=0+ 



(17) 



The transmissivity of a spin-cr electron through the mul- 
tiferric tunnel junctions reads 



K 



K 



(18) 



For a small applied bias voltages, the conductance G is 
determined by the states at the Fermi energy Ef i^^iSjie ^ 



d k|| dOr 
h J (27r)2 2tt 



■T^{EF,'k 



(19) 



FiglUJtop) shows the TAMR angular dependence ~ 
(1 - cos2(/.), at Ef = 5.5eV, A = 2eV, Vq = O.beV, 
and d = 2nm. It is clear that TAMR has C2v symme- 
try. For small J the spin-orbit interaction dominates the 
tunneling properties, we have positive TAMR. As J in- 
creases, the size of TAMR is influenced by the interplay 
between the exchange field and the induced spin-orbit 
interaction. The TAMR is typically ~ 0.1%, which is 
on the same order as in the Fe/GaAs/Au tunnel junc- 
tions have been recently realized experimentally^'^^. A 
transition from positive to negative TAMR is observed, 



which is consistent with the previous finding, i.e. Eq. ljlOp 
within the phenomenological model. On the other hand, 
the helicity of the spiral magnetic order q in helimagnetic 
multiferroics is experimentally controllable by a small 
(~ IkV/cm) transverse electric fieldi^. Consequently, 
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FIG. 2. The angular-dependence of TAMR on different 
strengths of the exchange field J{eV) and on the coefficient 7 
that enters the TAMR due to the competition between the ex- 
change field and the induced spin-orbit interaction. The used 
numerical values are Ef = 5.5cl/, A = 2eV , Vo = O.heV , and 
d = 2nm. Top: q = 0.28. Bottom: J = 0.1eV. 



we may electrically tune the spin-orbit coupling strength 
and thus TAMR in the helimagnet tunnel junctions (see 
Figl2;bottom)). 

Conclusions.- We studied the electron tunneling prop- 
erties through helical multiferroic junctions. The spiral 
magnetic ordering and the induced spin-orbit interaction 
in the multiferroic barrier lead to the TAMR effect. 
Due to the magnetoelectric coupling, the strength of the 
induced spin-orbit coupling is electrically controllable 
which reders a tunable TAMR by an external electric 
field. 

This research is supported by the DFG (Germany) un- 
der SFB 762. 
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